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Dark matter in the form of weakly interacting massive particles is predicted to become gravita- 
tionally captured and accumulate in stars. While the subsequent annihilations of such particles lead 
to the injection of energy into stellar cores, elastically scattering dark matter particles do not gener- 
ally yield enough energy to observably impact stellar phenomenology. Dark matter particles which 
scatter inelastically with nuclei (such that they reconcile the annual modulation reported by DAMA 
with the null results of CDMS and other experiments) , however, can be captured by and annihilate 
in compact stars at a much higher rate. As a result, old white dwarf stars residing in high dark 
matter density environments can be prevented from cooling below several thousand degrees Kelvin. 
Observations of old, cool white dwarfs in dwarf spheroidal galaxies, or in the inner kiloparsec of the 
Milky Way, can thus potentially provide a valuable test of the inelastic dark matter hypothesis. 
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I. INTRODUCTION 

Over the past two years, the DAMA/LIBRA collabo- 
ration has repeated their claim to have observed an an- 
nual modulation in their event rate, and have interpreted 
this as evidence of dark matter particles interacting with 
their detector at a confidence level of 8.2a fl . Such an 
annual modulation in the rate of dark matter scattering 
events is predicted to arise as a result of the variation 
in the velocity of the Earth with respect to the Milky 
Way's dark matter halo. Both the reported amplitude 
and phase of DAMA's signal are consistent with a dark 
matter interpretation. 

In typical dark matter models capable of producing the 
signal reported by DAMA, however, one would have ex- 
pected corresponding signals to appear in other direct de- 
tection experiments. To the contrary, several such exper- 
iments, including CDMS . CBESSTjJ, CoGeNT [1, 
XENON |L COUPP [|, ZEPLIN 07 Edelweiss 0, 
and KIMS [9J, have reported null results, strongly con- 
straining the dark matter interpretation of DAMA's an- 
nual modulation [l(|. This tension between DAMA and 
other direct detection experiments has motivated a num- 
ber of scenarios in which dark matter particles scat- 
ter preferentially with the sodium iodide (Nal) detec- 
tors of DAMA/LIBRA relative to the materials used in 
other experiments (such as germanium and silicon used 
in CDMS). Such scenarios feature dark matter particles 
which interact with nuclei through a resonance (Tlj . in- 
teract with nuclei with a momentum dependence causing 
them to scatter more efficiently with Nal than other tar- 
gets [IH, or which interact with nuclei largely through 
inelastic processes (l3| . In particular, dark matter par- 
ticles which, instead of elastically scattering with nuclei, 
scatter only through processes in which they are excited 
to a slightly (^100 keV) heavier state have generated a 
great deal of interest, and have been shown to be capable 
of reconciling the claims of the DAMA collaboration with 
the null results of other direct detection experiments [l3[ . 



In this article, we consider inelastic dark matter and 
study the role that such particle would play in compact 
stars, such as white dwarfs. In the case of dark mat- 
ter candidates which scatter elastically with nuclei, their 
capture by and annihilation in stars is expected to pro- 
vide an observable quantity of energy only in the most 
extreme environments, such as in a density spike of dark 
matter in the proximity of the supermassive black hole at 
the center of the Milky Way [14], or in the early phases 
of evolution of very massive population III stars [15| . In- 
elastic dark matter, in contrast, can potentially scatter 
with nuclei far more efficiently, but only if in posses- 
sion of sufficient momentum to overcome the mass split- 
ting between the ground and excited states. Because 
of the depth of the gravitational wells surrounding white 
dwarfs, inelastic dark matter particles will be accelerated 
to high velocities as they fall toward such a star, allowing 
them to easily accumulate enough momentum to scatter 
inelastically with the carbon or oxygen nuclei which make 
up the majority of the star's mass. As a result, inelas- 
tic dark matter particles will scatter in and be captured 
by white dwarf stars at very high rates. In their subse- 
quent annihilations, inelastic dark matter particles can 
deposit a significant amount of energy into these stars, 
potentially increasing their luminosity and/or tempera- 
ture. In particular, old and cool white dwarf stars that 
reside in regions of high dark matter density (such as in 
dwarf spheroidal galaxies, or in the inner regions of the 
Milky Way) may provide a valuable test of the inelastic 
dark matter scenario. 



II. CAPTURE OF INELASTIC DARK MATTER 
IN WHITE DWARF STARS 

More than two decades ago, the rate at which elas- 
tically scattering dark matter particles are captured in 
the Sun was calculated 16]. In many respects, the cap- 
ture rate of inelastically scattering dark matter particles 
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is similar. The key difference, however, is that in the in- 
elastic scenario the dark matter-nucleus scattering cross 
section depends on the velocity of the incoming particle. 
Above the threshold for inelastic scattering, the cross sec- 
tion is given by 



^inelastic — ^clastic A / 1 ni (1) 

V P v 

where S is the mass splitting between the dark matter 
and its excited state, /i is the dark matter-nucleus re- 
duced mass, and v is dark matter particle's velocity. The 
mass splitting, S, plays an important role in both direct 
detection experiments, and in the capture of dark matter 
particles onto the Sun. In particular, dark matter parti- 
cles with a splitting of 6 ~ 100 keV and a mass of 100 
GeV must have a velocity of v > 575 km/s if they are to 
inelastically scatter off of an iodine nucleus (A = 127). 
Including the motion of the Solar System through the 
dark matter halo, a non- negligible fraction of the dark 
matter particles are expected to exceed this velocity, es- 
pecially in the summer months, leading to the prediction 
of a high event rate with strong annual modulation in 
an experiment such as DAMA. In an experiment such 
as CDMS, which uses Germanium (A — 73) and Silicon 
(A = 28) target nuclei, however, the velocity threshold is 
significantly higher (v > 665 km/s), leading to consider- 
ably suppressed event rates. This relative suppression of 
the scattering rate with light target nuclei is the essence 
of the motivation for inelastic dark matter. 

In order for inelastic dark matter to avoid exceeding 
the limits of CDMS and XENON10, while also generating 
the annual modulation reported by DAMA, the scatter- 
ings must excite the dark matter into a state S ~ 80 — 150 
keV more massive that the ground state 0, [l7|- As a 
result of this mass splitting, inelastic dark matter scat- 
tering of off nuclei in the Sun is highly suppressed. In 
particular, essentially no scattering occurs in the Sun be- 
tween inelastic dark matter and nuclei lighter than car- 
bon. Although iron nuclei constitute only 0.16% of the 
Sun's mass, they dominate the capture rate of inelastic 
dark matter in the Sun. 

This suppression is not found, however, in the case of 
more dense stars, such as white dwarfs or neutron stars. 
The velocity of a dark matter particle that falls into a 
white dwarf from rest will exceed 5000 km/s at the sur- 
face, and 8000 km/s near the core (in contrast to approx- 
imately 1300 km/s near the core of the Sun). In this sit- 
uation, 28 //iv 2 <C 1, and the effects of the dark matter's 
mass splitting are negligible. For example, an inelastic 
dark matter particle falling from rest into a white dwarf 
consisting mostly of carbon or oxygen will inelastically 
scatter with a cross section within a few percent of the 
value found in the elastic limit. 

To calculate the capture rate of inelastic dark matter 
particles in a white dwarf, we follow the approach de- 
scribed in Refs. [H,[l!| for the case of capture in the Sun. 
For the Sun, we find results which are in good agreement 



with these previous studies (see the dashed lines in the 
right frame of Fig. [1}. 

White dwarf stars can be well modeled as a system 
in hydrostatic equilibrium with an equation of state of a 
cold fermi gas. Following Refs. [1(J[22], we describe the 
equation of state of a white dwarf as 

P = Cxf(x), 
P = C 2 x 3 , 

x = kf/(m e c), (2) 

where P is the pressure, p is the density of the star, and 
kf is the Fermi momentum. The numerical coefficients 
and the function f{x) are given by 

4 5 

8nm 3 e c 3 

f(x) = x(2x 2 -3){x 2 + 1) 1 ' 2 + Z\v{x+{l + x 2 ) 1/2 ],{Z) 

where h is the Planck constant. Substituting into the 
above equations for pressure and density, the equation of 
hydrostatic equilibrium becomes 

Ci l_d_ {r 2 df{x)\ = _ 4nGC2x3 (4) 
C2 r 2 dr \x 3 dr J 

which we solve numerically using a publicly available 
code 12211 that also includes corrections due to electro- 
static [2l| and general relativistic effects. This approach 
gives an accurate description of the stellar structure of 
white dwarfs and fits the data exquisitely (see, for exam- 
ple, Fig. 3 of Ref. [12). 

One way in which the calculation of inelastic dark mat- 
ter capture by white dwarfs departs from that described 
m Refs. [H, El is in the treatment of the star's opti- 
cal depth. In particular, for cross sections as large as 
are required to produce the DAMA annual modulation 
(cr~ 10~ 40 cm 2 ), dark matter particles are generally ex- 
pected to scatter many (typically several tens of) times 
as they pass through a white dwarf. In this optically 
thick limit, the rate at which dark matter particles are 
captured simplifies considerably to [23|, [24[ 

f 8n\ 3Gi?WD J^WD PDM /_% 

r c ~ -5- = , (5) 

where i?wH and M\vd are the radius and mass of the 
white dwarf, respectively, pdm is the density of dark mat- 
ter in the star's environment, and v is the dark matter's 
velocity dispersion (prior to infall into the gravitational 
potential of the star). Throughout our study, we will 
adopt M W d = 0.7M Q and i? W D = 0.0093i? Q as our 
default values (the mass distribution of white dwarfs is 
strongly peaked near 0.5-0.7 M Q [25j). Note that as long 
as we remain within the optically thick limit, the capture 
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FIG. 1: The rate at which energy is injected into white dwarf stars (solid) and into the Sun (dashed) by annihilating inelastic 
dark matter, as a function of the environment's dark matter density. In the left frame, we show results for a dark matter 
velocity dispersion of 20 km/s (appropriate for globular clusters [2{|), whereas in the right frame we use a velocity dispersion 
of 250 km/s and a star's velocity of 220 km/s (approprite for the Solar System). In each frame, upper (black) and lower (red) 
curves use dark matter- nucleon scattering cross sections of 10 -40 cm 2 and 10 -44 cm 2 , respectively (a ~ 10 -40 cm 2 is required 
to generate the DAMA annual modulation in the inelastic dark matter scenario, whereas a ~ 10 -44 cm 2 is approximately the 
maximum value allowed for elastically scattering dark matter candidates). In each case, we have adopted a dark matter mass 
of 100 GeV, and a mass splitting of S = 100 keV, although our results are not strongly sensitive to these values. The horizontal 
dotted lines denote the energy deposition rate above which observations of cool (3000-4000 K) white dwarfs would conflict with 
the predictions of inelastic dark matter. 
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rate does not depend on the interaction cross section of 
the dark matter, but only on its number density and ve- 
locity distribution. 

The large cross section required in the inelastic dark 
matter scenario quickly leads to an equilibrium between 
the rates of capture and annihilation in a white dwarf, as- 
suming that the dark matter particles are capable of self- 
annihilation (which might not be true if the dark matter 
carries a conserved quantum number rather being stabi- 
lized by a parity, for example). As a result, the capture 
rate can be converted into the rate at which additional 
energy is injected into the star's core. In a typical white 
dwarf star (i? WD « 0.0093i? Q and M WD « O.7M ), in- 
elastic dark matter annihilations provide a contribution 
to the star's luminosity that is given by: 

L PS r c 771DM 

« 2.5 x 10 28 GeV/s ( PUM , ) ( 2 J^ll) 
V 1 GeV/cm 3 /V « / 

« 4xlO*Wsf PDM 3 V 220km/S V (6) 
VlGcV/cm 3 / V v J W 

In Fig. [I] we show the rate at which inelastic dark mat- 
ter deposits energy into a white dwarf (or into the Sun) , 
as a function of the environment's dark matter density. 
We show results for dark matter with velocity dispersions 
of 20 km/s (left) and 250 km/s (right), and use two val- 
ues for the dark matter-nucleon cross section: 10 -44 cm 2 
(lower red curves) and 10 -40 cm 2 (upper black curves). 
In the case of the larger (smaller) cross section choice, a 



white dwarf is within the optically thick (thin) limit. As 
a cross section of a « 10 -40 cm 2 is required to generate 
the DAMA modulation, we will focus on the optically 
thick case in the remainder of the study. 

The coldest and least luminous observed white dwarfs 
have temperatures as low as 3000 K and luminosities of 
(2 — 3) x 10 28 erg/s (following from the Stefan-Boltzmann 
law, a 3000 K blackbody of radius 0.01 Rq is predicted 
to have a luminosity of 2.8 x 10 28 erg/s). From these ob- 
servations, we can conclude that any additional source of 
energy (such as dark matter annihilations) must not ex- 
ceed approximately 2 x 10 28 erg/s. Within the context of 
inelastic dark matter, this translates to an approximate 
limit of pdm £ 500 GeV/cm 3 (p D M £ 40 GeV/cm 3 ) for 
v = 250 km/s (v = 20 km/s). Although these density 
limits are orders of magnitude higher than the average 
dark matter density in the local reg ion of the Milky Way 
(pdm, local ~ 0.3-0.5 GeV/cm 3 [27|), there exist regions 
that are thought to contain sufficiently high quantities of 
dark matter (and/or with sufficiently slow velocity dis- 
tributions) to make white dwarf stars useful probes of 
inelastic dark matter. We discuss some of these possibil- 
ities in the following section. 



III. DARK MATTER DENSITIES IN REGIONS 
CONTAINING WHITE DWARF STARS 

Many of the oldest (and coldest) known white dwarf 
stars are located in the dense collections of stars known 
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as globular clusters. If high densities of dark matter 
were also present in such systems, these very old and 
cool white dwarfs would provide an ideal environment to 
search for the effects of inelastic dark matter. Despite the 
very high densities of baryonic matter in globular clusters 
(often containing hundreds of thousands of stars within 
a volume smaller than a cubic parsec), however, there 
is no compelling dynamical evidence that these systems 
contain significant quantities of dark matter. 

Within the modern paradigm of globular cluster for- 
mation, these objects form in the disk of galaxies rather 
than in the center of dark matter halos (in contrast to 
dwarf spheroidal galaxies, for example). Indeed, observa- 
tions of the Antennae Galaxy dramatically demonstrate 
the formation of globular clusters within its disk. Re- 
gardless, the formation of a globular cluster is generally 
expected to enhance the dark matter density through 
adiabatic contraction. In particular, the orbits of dark 
matter particles that are gravitationally bound to the 
cluster become contracted as the baryons (ie. the stars) 
themselves gather more densely. This can lead to a con- 
tracted dark matter density in the globular cluster's core 
approximately given by 



PDM,c 



Pb,c 



Pdm, 

Pb,i 



( Pism \ 

\ Pcorc / 



1/4 



(7) 



where pb.corc is the density of baryons in the globular 
cluster's core, pDM.i is the initial density of gravitation- 
ally bound dark matter particles prior to adiabatic con- 
traction, is the initial density of baryons prior to 
adiabatic contraction, pism is the total density of the in- 
terstellar medium, and p C orc is the total density in the 
core of the cluster. This result was found analytically 
using the standard adiabatic contraction model of Blu- 
menthal et al. 28], and was verified numerically using 
the adiabatic contruction code, Contra (29j . 

Considering as an example the nearby globular cluster 
M4, we would expect that within an initial distribution of 
dark matter with a velocity dispersion of v ~ 220 km/s, 
only ~ 10~ 3 — 10~ 4 or less of the particles will be grav- 
itationally bound, thus significantly limiting the impact 
of adiabatic contraction. After adiabatic contraction, we 
estimate that the dark matter density in the core of M4 
could be enhanced to on the order of a few GeV/cm 3 
which, although an order of magnitude greater than prior 
to the formation of the cluster, is well below the density 
needed to significantly effect white dwarf stars. 

As an alternative to globular clusters, we also con- 
sider white dwarfs within the dark matter halos of dwarf 
spheroidal galaxies. Dwarf spheroidals are highly dense, 
dark matter dominated systems, and contain very old 
stars (30j . making them excellent environments in which 
to study the impact of inelastic dark matter on white 
dwarfs. Ref. [31) . for example, describes seven dwarf 
spheroidals (Carina, Draco, Fornax, Leo I, Leo II, Sculp- 
tor, and Sextans) which have dark matter densities of 
40-150 GeV/cm 3 and velocity dispersions of ^10-20 km/s 
within their inner 10-20 parsecs. If old and cold white 



dwarf stars are present in such environments, we pre- 
dict that inelastic dark matter will inject energy into 
them at a rate rivaling their luminosity. As a result, 
the observation of a white dwarf in such an environment 
with a temperature below approximately 4000 K would 
conflict with the inelastic dark matter hypothesis. Al- 
though dwarf spheroidals are typically very distant, and 
thus challenging to observe, there are 3-4 known dwarf 
spheroidals (Segue I, Segue II, Ursa Major II, and pos- 
sibly Wilman I) within ^30 kpc of the Solar System, 
making them within the approximate range of a 30 me- 
ter telescope [32J. In the more distant future, a 100 meter 
telescope could place white dwarf stars in all Milky Way 
dwarf spheroidals within reach. 

Lastly, the dark matter density is generally predicted 
to increase as one approaches the center of the Milky 
Way. According to the frequently used Navarro-Frenk- 
White (NFW) profile (33[, for example, the smooth com- 
ponent of the dark matter density in the Milky Way can 
be parameterized by 



1 



Pdm oc 



r[l + {r/R s )f 



(8) 



where r is the distance to the Galactic Center, and 
R s s=a 20 kpc is the scale radius. In addition, the ve- 
locity dispersion is predicted to decrease according to 
{r/Rs) 1 ' 2 . Comb ining these effects, we estimate 
that the rate at which inelastic dark matter will accu- 
mulate in a white dwarf located at 1 kpc (100 pc) from 
the Galactic Center is ~40 (~1250) times greater than 
the rate predicted in the vicinity of the Solar System. 
Adopting an NFW profile, we expect a sizable fraction 
of the luminosity of cold white dwarf stars within the in- 
ner ^100 pc of the Milky Way to originate from inelastic 
dark matter annihilations. This effect can be consider- 
ably enhanced if the Milky Way's dark matter halo were 
adiabatically contracted by baryons 34|. We therefore 
expect future infrared surveys of the inner Milky Way to 
provide a valuable test of inelastic dark matter. 



IV. CONCLUSIONS 

In this paper, we have studied the capture and sub- 
sequent annihilation of inelastic dark matter particles in 
white dwarf stars. We find that inelastic dark matter (as 
motivated to generate the annual modulation claimed by 
DAMA/LIBRA without exceeding the limits from other 
direct detection experiments) is predicted to be captured 
by compact stars at a rate multiple orders of magni- 
tude higher than elastically scattering dark matter. This 
makes inelastic dark matter a very efficient source of en- 
ergy for white dwarf stars, potentially leading to observ- 
able consequences. 

The role of inelastic dark matter is most pronounced 
in those stars that reside in regions with a high density 
(and/or low velocity dispersion) of dark matter. In par- 
ticular, we find that very old white dwarf stars in the 
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inner tens of parsecs of dwarf spheroidal galaxies will be 
powered primarily by inelastic dark matter annihilations, 
preventing their otherwise expected cooling. This is also 
predicted to be the case for old white dwarfs found within 
the inner kiloparsec of the Milky Way. If future surveys 
were to discover very cold (T < 4000 K) white dwarfs in 
environments such as these, it would disfavor the inelas- 
tic dark matter hypothesis. In contrast, the lack of such 
stars discovered in future observations could be used to 
bolster the case for inelastic dark matter. 
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